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The mass of the 1 _+ exotic meson, created with hybrid interpolating fields, is explored in numer- 
ical simulations of quenched QCD on large ( 20 3 x 40 ) lattices to obtain good control of statistical 
and finite volume errors. Using the Fat-Link Irrelevant Clover (FLIC) fermion action, the proper- 
ties of the 1 v are investigated at light quark masses approaching 25 MeV (m n /m p ~ 1/3). Under 
the standard assumption that the coupling to the quenched a\v[ channel comes with a negative 
metric, our results indicate that the 1 h exotic exhibits significant curvature close to the chiral 
limit, suggesting previous linear extrapolations have overestimated the mass of the 1" + . We find 
for the first time in lattice studies a 1 mass consistent with the 7Ti(1600) candidate. We also find 
a strangeness ±1 J p = 1~ state with a mass close to 2 GeV. 
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I. INTRODUCTION 

The masses of the so-called 'exotic' mesons are attract- 
ing considerable attention from the experimental commu- 
nity [EHIESHISSIlasa vehicle for the elucidation of 
the relatively unexplored role of gluons in QCD. The Par- 
ticle Data Group |£| reports two candidates for the 1 h 
exotic, the tti(1400) at 1.376(17) GeV, and the 7Ti(1600) 
at 1.596^4 GeV. The experimental status of these states 
is an issue that continues to attract attention jy, |7j . 

Lattice QCD provides a first principles approach to 
nonperturbative calculations of QCD, indispensable in 
determining the hadron mass spectrum. In the early 
work of Ref. ^(j, the UKQCD Collaboration made use 
of gauge-invariant non-local operators to explore P and 
D-wave mesons, as well as exotics. They used a tadpole- 
improved clover action, with 375 configurations on a 
16 3 x 48 lattice. Their calculation was performed at a 
single quark mass corresponding to approximately that 
of the strange quark and reported a 1 h exotic mass of 
1.9(4) GeV. 

In 1997, the MILC Collaboration used local operators 
formed by combining the gluon field strength tensor and 
standard quark bilinears |ll| , the same approach we take 
in this paper. The highly anisotropic lattices employed 
allowed many time slices to be used in determining the 
mass of the exotic. 20 3 x 48 and 32 3 x 64 lattices with mul- 
tiple fermion sources per lattice were considered. Using 
the Wilson action, they suggested a 1 h mass of 1.97(9) 
GeV from a linear extrapolation in 1/k from m 2 > 0.64 
GeV. 

The SESAM collaboration analyzed dynamical 
fermion configurations at four quark masses correspond- 
ing to rrv^ > 0.455(9) GeV 2 . Their linear extrapolation 
resulted in 1.9(2) GeV for the mass of the 1 h exotic 
where the scale was determined via the Sommer param- 
eter [13 r via Ref. [l|. 



Further work using the Clover action, but this time 
with local interpolators was performed by Mei et al. in 
Ref. 0] . Very heavy quark masses were used to get good 
control of statistical errors. Their linear extrapolation 
from to 2 > 1.05 GeV 2 suggested a mass of 2.01(7) GeV. 

In 2002 the MILC Collaboration published new work 
|l6j using dynamical improved Kogut-Susskind fermions 
on 20 3 x 48 and 32 3 x 64 lattices. Linearly extrapolating 
from m 2 > 0.488 GeV 2 , Bernard et al. quote two sets of 
results for the 1 h mass. Using r\ = 0.34 fm and ^[g 
[■[() MeV they report a 1~+ mass of 1.85(7) and 2.03(7) 
GeV respectively. 

Michael 01 provides a good summary of work to 2003, 
concluding that the light-quark exotic is predicted by lat- 
tice studies to have a mass of 1-9(2) GeV. 

In order to minimize the need for extrapolation one 
requires access to quark masses near the chiral regime 
on large physical volumes. Our study considers a phys- 
ical volume of (2.6 f m) 3 , and the 0(a)-improved FLIC 
fermion action [LSI llii, |20| whose improved chiral proper- 
ties pH permit the use of very light quark masses which 
are key to our results. 



II. VECTOR MESONS ON THE LATTICE 

Consider the momentum-space meson two-point func- 
tion for t > 0, 

G%(t,p) = Y^e- iff - S (n\ti(t,£)xi\0,0m) (1) 

X 

where i,j label the different interpolating fields and /i, v 
label the Lorentz indices. At the hadronic level, 
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TABLE I: J quantum numbers and their associated meson interpolating fields. 
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where the \E,p' , s) are a complete set of hadronic states, 
of energy E, momentum p', and spin s, 

\E,p', 3 )(E,p',s\=I . (2) 

E,j? ,s 

We denote the vacuum couplings as follows: 

(n\^\E,p',s) = \ E e^p\s) 

(E,p' 7 s\xC\n) = \ j E *el{p\s) , 

where the four- vector p ' = (E,p') is introduced. 

We can translate the sink operator from a; to to write 
this as 

£ e-^fll 4(0) e iPs - At \E,p', s) 
x (£,^',51x^(0)10) 
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A l e M (p, s)X J e*(p,s) 



(3) 



We shall label the states which have the \ interpolating 
field quantum numbers as \a) for a — 1, 2, • • ■ , AT, thus 
replacing ^2 E with £^ a . In general the number of states, 
N, in this tower of excited states may be very large, but 
we will only ever need to consider a finite set of the lowest 
energy states here, as higher states will be exponentially 
suppressed as we evolve to large Euclidean time. Finally, 
the transversality condition: 



E 



(4) 



implies that for p = 0, we have 
<4(i,0)=0 

G%(t,d)=Y,5kiK\Ce~ m « t ■ (5) 

a 

Since G^, G^and G 3 g are all estimates for the same 
quantity we add them together to reduce variance, form- 
ing the sum 

siij _ fiij , fiij , pi] 



Evolving to large Euclidean time will suppress higher 
mass states exponentially with respect to the lowest-lying 
state, leading to the following definition of the effective 
mass 



M«(t)= In 



G ij (t,0) 
Gv(t + 1,0) 



(6) 



The presence of a plateau in M e ff as a function of time, 
then, signals that only the ground state signal remains. 



III. LATTICE SIMULATIONS 
A. Interpolating Fields 

The formulation of effective interpolating fields for the 
creation and annihilation of exotic meson states contin- 
ues to be an active area of research. For example, one 
can generalize the structure of the interpolating fields to 
include non-local components where link paths are in- 
corporated to maintain gauge invariance and carry the 
nontrivial quantum numbers of the gluon fields |l0l Il2| . 
In this case, numerous quark propagators are required 
for each gauge field configuration rendering the approach 
computationally expensive. The use of non-local interpo- 
lating fields affords greater freedom in creation of opera- 
tors and facilitates access to exited states through varia- 
tional techniques, but does not lead to an increase in sig- 
nal for the ground state 1 ^ exotic commensurate with 
the increased computational cost of this approach. Ex- 
otic quantum numbers may also be obtained from four- 
quark (qqqq) operators, but in practice, these tend to 
exhibit larger statistical fluctuations pTol fill fl2l| . 

We consider the local interpolating fields summarized 
in Table [I] Gauge-invariant Gaussian smearing [2^, |2^| 
is applied at the fermion source (t = 8), and local sinks 
are used to maintain strong signal in the two-point cor- 
relation functions. Chromo-electric and -magnetic fields 
are created from 3-D APE-smeared links [24J at both 
the source and sink using the highly- improved 0(a 4 )- 
improved lattice field strength tensor |25| described in 
greater detail below. 

While all four 1 h interpolating fields of Table [I] have 
been used to create the 1 h state on the lattice, in this 
work we will focus on the results for the two 1 ^ inter- 
polating fields coupling large spinor components to large 
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spinor components: 

X2 
X3 



ie jk iq a lk Bfq b 
ie ]k iq a lilkBt b q b . 



(7) 



These interpolating fields provide the strongest signal for 
the 1 h state. In both fields the gluonic component is 
contributing axial-vector quantum numbers. 



B. Lattice Field Strength Tensor 

In order to obtain the chromo-electric and chromo- 
magnetic fields with which we build the hybrid operators, 
we make use of a modified version of APE smearing, in 
which the smeared links do not involve averages which 
include links in the temporal direction. In this way we 
preserve the notion of a Euclidean 'time' and avoid over- 
lap of the creation and annihilation operators. 

In this study, the smearing fraction a = 0.7 (keeping 
0.3 of the original link) and the process of smearing and 
ST/ (3) link projection is iterated six times [2(|. Each iter- 
ation of our modified APE-smearing algorithm proceeds 
as 

Ui(x) (l-a)E7i(a!) 

3 



3 = 1 
3 



4 - ) U H X - 3) U & ~ 3') Uj(x - j + I) 

3=1 



U 4 {x) -> (l-a)U 4 {x) 



a 



3=1 

3 



U] (x - j) U 4 (x - j) Uj (x - ] + 4) 



3=1 



Smearing the links permits the use of highly improved 
definitions of the lattice field strength tensor, from which 
our hybrid operators are derived. Details of the C(a 4 )- 
improved tensor are given in |25j . This amount of smear- 
ing is suitable for the creation of exotic mesons. 
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FIG. 1: Exotic meson propagator for interpolator X2- Results 
are shown for every 2nd quark mass in the simulation. Lower 
lines correspond to heavier quark masses. For all but the 
heaviest mass, the signal is lost after t=12. 



-4 - 



-6 - 



10 



■12 




FIG. 2: Exotic meson propagator for interpolator X3- Results 
are shown for every 2nd quark mass in the simulation. Lower 
lines correspond to heavier quark masses. 



with clover actions j2J , and minimizes the effect of renor- 
malization on the action improvement terms |27| . Access 
to the light quark mass regime is enabled by the improved 
chiral properties of the lattice fermion action [2l|. By 
smearing only the irrelevant, higher dimensional terms in 
the action, and leaving the relevant dimension-four oper- 
ators untouched, short distance quark and gluon interac- 
tions are retained. Details of this approach may be found 
in reference [Tsj . FLIC fermions provide a new form of 
nonperturbative O(a) improvement |2lL |27| where near- 
continuum results are obtained at finite lattice spacing. 



Propagators are generated using the fat-link irrelevant 
clover (FLIC) fermion action where the irrelevant 
Wilson and clover terms of the fermion action are con- 
structed using fat links, while the relevant operators use 
the untouched (thin) gauge links. Fat links are created 
via APE smearing |2J|. In the FLIC action, this reduces 
the problem of exceptional configurations encountered 



D. Gauge Action 

We use quenched-QCD gauge fields created by the 
CSSM Lattice Collaboration with the 0(a 2 ) mean-field 
improved Liischer-Weisz plaquette plus rectangle gauge 
action 28] using the plaquette measure for the mean link. 
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FIG. 3: Effective mass for interpolator \2- Plot symbols are 
as for the corresponding propagator plot. 
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FIG. 4: As for Fig. but for interpolator X3- Signal is lost 
after t — 11. 



The gauge-field parameters are defined by 
S e = b 4 E^ReTta-iVCr)) 

V>fl 

- ^ X! ^Tr(2-^(x)) , 

x fii/ 
v>H 

where P M „ and i? M „ are defined in the usual manner and 
the link product i? M „ contains the sum of the rectangular 
1x2 and 2x1 Wilson loops. 

The CSSM configurations are generated using the 
Cabibbo-Marinari pseudo-heat-bath algorithm [2^| using 
a parallel algorithm with appropriate link partitioning 
|30| . To improve the ergodicity of the Markov chain pro- 
cess, the three diagonal SU(2) subgroups of SU(3) are 
looped over twice |3l| and a parity transformation [32| is 
applied randomly to each gauge field configuration saved 
during the Markov chain process. 



FIG. 5: Effective mass for the pseudoscalar pion interpolator 
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FIG. 6: Effective mass for the axial-vector pion interpolator 
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E. Simulation Parameters 

The calculations of meson masses are performed on 
20 3 x 40 lattices at (3 = 4.53, which provides a lattice 
spacing of a — 0.128(2) fm set by the Sommer parameter 
r = 0.49 fm. 

A fixed boundary condition in the time direction is 
used for the fermions by setting Ut(£, N t ) = V x in 
the hopping terms of the fermion action, with periodic 
boundary conditions imposed in the spatial directions. 

Eight quark masses are considered in the calculations 
and the strange quark mass is taken to be the third heavi- 
est quark mass. This provides a pseudoscalar mass of 697 
MeV which compares well with the experimental value of 
(2M| - Ml) 1 / 2 = 693 MeV motivated by leading order 
chiral perturbation theory. 

The analysis is based on a sample of 345 configura- 
tions, and the error analysis is performed by a third-order 
single-elimination jackknife, with the y 2 per degree of 
freedom (x 2 /(io/) obtained via covariance matrix fits. 
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FIG. 7: Effective mass for the interpolator \2 with a strange 
quark. 
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FIG. 8: As for Fig.[7| but for interpolator X3- 

IV. RESULTS 

Figures ^ and El show the natural log of the correla- 
tion functions calculated with interpolators %2 and X3 
from Eq. [7] respectively. The curves become linear af- 
ter two time slices from the source, corresponding to ap- 
proximately 0.256 fm. This is consistent with Ref. . 
where a similar effect is seen after approximately 3 to 4 
time slices, corresponding to 0.21 to 0.28 fm following the 
source. 

Figures |21 and 0] show the effective mass for the two 
different interpolators. For clarity, we have plotted the 
results for every second quark mass used in our simula- 
tion. The plateaus demonstrate that we do indeed see 
an exotic signal in quenched lattice QCD. This is signifi- 
cant, as we expect the two interpolating fields to possess 
considerably different excited-state contributions, based 
on experience with pseudoscalar interpolators |33| . 

For example, the approach to the pion mass plateau is 
from above (below) for the pseudo-scalar (axial-vector) 
interpolating field as illustrated in Fig. [5] (Fig. . This 
exhibits the very different overlap of the interpolators 



FIG. 9: A survey of results in this field. The MILC results 
are taken from fill] and show their Q 4 , 1 _+ — » U + results, 
fitted from t = 3 to t = 11. Open and closed symbols denote 
dynamical and quenched simulations respectively. 

with excited states. As in the 1 h interpolators, the role 
of 74 in the pion interpolators is to change the sign with 
which the large-large and small-small spinor components 
are combined. We also present results for the strangeness 
±1 analogue of the 1 h in Figs. Hand |B1 

TablelTTlsummarizes our results for the mass of the 1 
meson, with the squared pion-mass provided as a measure 
of the input quark mass. The agreement observed in 
the results obtained from the two different 1 h hybrid 
interpolators provides evidence that a genuine ground- 
state signal for the exotic has been observed. 

Table IIIII summarizes our results for the mass of the 
strangeness ± 1, J p = 1~ meson. 

Finally, in Fig. [|]] we summarize a collection of results 
for the mass of 1 h obtained in lattice QCD simulations 
thus far. The current results presented herein are com- 
pared with results from the MILC [H and SESAM 
|12| collaborations, both of which provide a consistent 
scale via ro. 

Our results compare favorably with earlier work at 
large quark masses. Agreement within one sigma is ob- 
served for all the quenched simulation results illustrated 
by filled symbols. It is interesting that the dynamical 
Wilson fermion results of the SESAM collaboration ^1 
tend to sit somewhat hi ghe r as this is a well known effect 
in baryon spectroscopy [l]| IM Li! • 



V. PHYSICAL PREDICTIONS 

In comparing the results of quenched QCD simulations 
with experiment, the most common practice is to simply 
extrapolate the results linearly in m q or to the phys- 
ical values. However, such an approach provides no op- 
portunity to account for the incorrect chiral nonanalytic 
behavior of quenched QCD El S H3 • 

Unfortunately, little is known about the chiral nonana- 
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TABLE II: 1 + Exotic Meson mass m (GeV) vs square of pion mass ml (GeV 2 



ml 




X 2 fit 10- 1 1 




X 2 fit 10-12 




X3 fit 10-11 




m 


X 2 /dof 


m 


X 2 /dof 


m 


X 2 /dof 



0.693(3) 
0.595(4) 
0.488(3) 
0.381(3) 
0.284(3) 
0.215(3) 
0.145(3) 
0.102(4) 



2.15(12) 
2.11(12) 
2.07(12) 
2.01(12) 
1.97(13) 
1.92(14) 
1.85(17) 
1.80(23) 



0.69 
0.77 
0.85 
0.91 
0.78 
0.78 
0.57 
0.13 



2.16(11) 
2.12(11) 
2.08(12) 
2.03(12) 
1.98(13) 
1.92(14) 
1.84(17) 
1.75(23) 



0.44 
0.51 
0.59 
0.65 
0.55 
0.40 
1.76 
3.04 



2.20(15) 
2.18(16) 
2.15(17) 
2.14(19) 
2.27(29) 
2.25(31) 
2.26(37) 
2.46(58) 



0.45 

0.46 

0.41 

0.29 

0.00012 

0.02 

0.02 

0.03 



TABLE III: Strangeness ±1 1 Meson mass m (GeV) vs square of pion mass m^(GeV 2 



ml 




X 2 fit 10-11 




X 2 fit 10-12 




X3 fit 10-11 




m 


X 2 /dof 


m 


X 2 /dof 


m 


X 2 /dof 



0.693(3) 
0.595(4) 
0.488(3) 
0.381(3) 
0.284(3) 
0.215(3) 
0.145(3) 
0.102(4) 



2.11(12) 
2.09(12) 
2.07(12) 
2.04(12) 
2.01(13) 
1.99(13) 
1.97(13) 
1.96(14) 



0.76 
0.81 
0.85 
0.88 
0.85 
0.87 
0.73 
0.56 



2.12(11) 
2.10(12) 
2.08(12) 
2.05(12) 
2.02(12) 
2.00(12) 
1.97(13) 
1.96(14) 



0.51 
0.55 
0.59 
0.63 
0.63 
0.64 
0.54 
0.39 



2.17(16) 
2.16(16) 
2.15(17) 
2.15(18) 
2.25(20) 
2.11(20) 
2.12(22) 
2.09(24) 



0.44 
0.44 
0.41 
0.36 
0.22 
0.29 
0.11 
0.01 



lytic behavior of the 1 h meson. Ref. 38] provides a full 
QCD exploration of the chiral curvature to be expected 
from transitions to nearby virtual states and channels 
which are open at physical quark masses. While vir- 
tual channels act to push the lower-lying single-particle 
1 h state down in mass, it is possible to have sufficient 
strength lying below the 1 h in the decay channels such 
that the 1 h mass is increased l39t JiOl . Depending on 
the parameters considered in Ref. [385 governing the cou- 
plings of the various channels, corrections due to chiral 
curvature are estimated at the order of +20 to —40 MeV. 

Generally speaking, chiral curvature is suppressed in 
the quenched approximation. For mesons, most of the 
physically relevant diagrams involve a sea-quark loop and 
are therefore absent [4(J,|4l]]. However, the light quenched 
rj meson can provide new nonanalytic behavior, with the 
lowest order contributions coming as a negative-metric 
contribution through the double-hairpin diagrams. Not 
only do these contributions alter the 1 h mass through 
self-energy contributions, but at sufficiently light quarks 
masses, open decay channels can dominate the two-point 
correlator and render its sign negative. 

For the quenched 1 v meson, the a\rj' channel can be 
open. Using the pion mass as the rj' mass a direct calcu- 
lation of the mass of an airj 1 two-particle state indicates 
that the 1 h hybrid lies lower than the two-particle state 
for heavy input quark mass. This indicates that the hy- 



brid interpolator is effective at isolating a single-particle 
bound state as opposed to the two-particle state at heavy 
quark masses. This is particularly true for the case here, 
where long Euclidean time evolution is difficult. 

As the light quark mass regime is approached, the 
trend of the one and two-particle states illustrated in 
Fig. EH suggests that they either merge or cross at our 
second lightest quark mass, such that the exotic 1 *~ 
may be a resonance at our lightest quark mass and at 
the physical quark masses. We note that the exotic 1 h 
mass displays the common resonance behavior of becom- 
ing bound at quark masses somewhat larger than the 
physical quark masses. This must happen at sufficiently 
heavy quark masses by quark counting rules, i.e 2q — -*■ 4g 
for the 1 h to air]' transition. 

One might have some concerns about airy' contamina- 
tions in the two-point correlation function affecting the 
extraction of the 1 ^ meson mass 01 • However we can 
already make some comments. 

Under the assumption that the coupling to the 
quenched air/' channel comes with a negative metric, as 
suggested by chiral perturbation theory arguments, and 
from the observation that our correlation functions are 
positive, then it would appear that our interpolators cou- 
ple weakly to the decay channel. Furthermore, at heavy 
quark masses the correlation function is dominated by 
the 1 h bound state already at early Euclidean times 
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FIG. 10: The 1 h exotic meson mass obtained from fits of the 
effective mass of the hybrid interpolator X2 from t — 10 — > 
12 (full triangles) are compared with the air]' two-particle 
state (open triangles). The extrapolation curves include a 
quadratic fit to all eight quark masses (dashed line) and a 
linear fit through the four lightest quark masses (solid line). 
The full square is result of linear extrapolation to the physical 
pion mass, while the open square (offset for clarity) indicates 
the 7ri(1600) experimental candidate. . 
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FIG. 11: Extrapolation of the associated strangeness ±1 
J p — 1~ state obtained from %2- Symbols are as in Fig. 1101 



suggesting that coupling to the decay channel is weak. 

Thus we conclude that the hybrid interpolating fields 
used to explore the 1 ^ quantum numbers are well-suited 
to isolating the single-particle 1 h exotic meson. 

Moreover, since the mass of the a\rj channel is simi- 
lar or greater than the single-particle 1 v state, one can 
conclude that the double-hairpin a\rj' contribution to the 
self energy of the single-particle 1 ^ exotic meson is re- 
pulsive in quenched QCD. Since the curvature observed 
in Fig. [5] reflects attractive interactions, we can also con- 
clude that quenched chiral artifacts are unlikely to be 
large. 

Hence we proceed with simple linear and quadratic ex- 
trapolations in quark mass to the physical pion mass, 
with the caution that chiral nonanalytic behavior could 



provide corrections to our simple extrapolations the order 
of 50 MeV in the 1 1 mass jig. 

Figures EH and ^2 illustrate the extrapolation of the 
1 h exotic and its associated strangeness ±1 1~ state 
to the limit of physical quark mass. We perform the 
linear fit using the four lightest quark masses and fit 
the quadratic form to all 8 masses. A third-order 
single-elimination jackknife error analysis yields masses 
of 1.74(24) and 1.74(25) GeV for the linear and quadratic 
fits, respectively. These results agree within one stan- 
dard deviation with the experimental 7Ti(1600) result of 
1.596±f| GeV, and exclude the mass of the 7Ti(1400) can- 
didate. 

The associated parameters of the fits are as follows. 
The linear form 

m!-+ = a + a 2 ml , 

yields best fit parameters of 

a a = 1.73 ±0.15 GeV, 
a 2 =0.85 ±0.35 GeV -1 . 



a 2 ml 



4 

04 777^ . 



PC 



V 



The quadratic fit, with formula 

777!-+ = a 

returns parameters 

no = +1.74 ±0.15 GeV, 
a 2 = +0.91 ±0.39 GeV" 1 
a 4 = -0.46 + 0.35 GeV" 3 

VI. CONCLUSION 

We have found a compelling signal for the J 
exotic meson , from which we can extrapolate a physi- 
cal mass of 1.74(24) GeV. Thus for the first time in lat- 
tice studies, we find a 1 ^ mass in agreement with the 
7Ti(1600) candidate. 

The X2 interpolating field appears to be extremely use- 
ful for avoiding contamination from the air] 1 channel, and 
thus is an excellent choice for this kind of study. 

We have also presented the first results for a 
strangeness ±1 partner of the exotic 1 h meson lying 
at 1.92(15) GeV. 

Looking forward, it will be important to quantify the 
effects of the quenched approximation. We plan to re- 
visit these calculations at some future point using full 
dynamical FLIC fermions [43l |44| . Of particular inter- 
est will be the extent to which the curvature observed in 
approaching the chiral regime is preserved in full QCD. 

Additionally, we intend to explore the dependence of 
the exotic signal on the nature of the fcrmion source. 
Whilst the rapidity with which we establish a plateau in 
our effective mass plots suggests that our current fermion 
operator smearing is near optimal for isolating the ground 
state, it might be possible to reduce the statistical errors 
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through a careful selection of parameters coming out of 
a systematic exploration of the parameter space. 

Finally, a detailed finite volume analysis should be per- 
formed in order to further explore the role of the two- 
body decay channel. 
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